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We show how to numerically calculate several quantities that characterize topological order starting from 
a microscopic fractional quantum Hall (FQH) Hamiltonian. To find the set of degenerate ground states, we 
employ the infinite density matrix renormalization group (iDMRG) method based on the matrix-product state 
(MPS) representation of FQH states on an infinite cylinder. To study localized quasiparticles of a chosen topo- 
logical charge, we use pairs of degenerate ground states as boundary conditions for the iDMRG. We then show 
that the wave function obtained on the infinite cylinder geometry can be adapted to a torus of arbitrary mod- 
ular parameter, which allows us to explicitly calculate the non-Abelian Berry connection associated with the 
modular T-transformation. As a result, the quantum dimensions, topological spins, quasiparticle charges, chiral 
central charge, and Hall viscosity of the phase can be obtained using data contained entirely in the entanglement 
spectrum of an infinite cylinder. 



Over the last decades several new kinds of phases have 
been discovered which cannot be characterized by sponta- 
neous symmetry breaking, but instead exhibit topological or- 
der (T). A prominent example of such a phase is the fractional 
quantum Hall (FQH) effect which occurs due to electron- 
electron interactions in a Landau level of fractional filling v 
0. Topologically ordered phases have unusual properties 
such as excitations which carry fractions of the elementary 
charge and have exotic exchange statistics {i.e., they are nei- 
ther fermions nor bosons) J3]|4]. Recent proposals that the 
FQH state at filling v = 5/2 |5| can be used to build a 'topo- 
logical' quantum computer have spurred both theoretical and 
experimental advances in FQH systems li6l n"Tl . 

Much of our current understanding of FQH phases results 
from trial wave functions, effective field theories, and numer- 
ical simulations. In particular, exact diagonalization (ED) is 
a powerful tool for probing the ground state and excitations 
of these phases |4][T2). Certain properties, such as the quan- 
tum dimensions of the quasiparticles, can then be extracted 
from just the ground-state wave functions |[T3l{T7l . However, 
the exponential growth of the Hilbert space as a function of 
system size means that ED is prohibitively expensive beyond 
w20 particles. While it is possible to obtain very accurate re- 
sults for some systems with a small correlation length, such as 
the Coulomb state at filling factor v = 1/3, more complicated 
systems (system with higher Landau level, hierarchy states, or 
non-Abelian phases) are much harder to access and finite size 
effects are much stronger. It would be advantageous to de- 
velop numerical techniques that allow for both larger system 
sizes and a more complete characterization of the topological 
order. 

In this paper we address the limitations of ED by using the 
infinite density matrix renormalization group (iDMRG) algo- 
rithm to obtain the matrix product state (MPS) representation 
of the ground-state of FQH Hamiltonians on infinitely long 
cylinders with finite circumference L (see Fig. [TJ. The space 
of MPS have been shown to be an exact and efficient repre- 
sentation of the model QH states (|T8l [19J; iDMRG extends 



these results to non-model states by variationally optimizing 
an iMPS with respect to a microscopic Hamiltonian ||20l . For 
a system of size L, the computational complexity reduces 
from 0(b L2 ) via ED to 0(b L ) for DMRG, with b > 1. Sev- 
eral groups have by now implemented finite system DMRG 
to simulate FQH Hamiltonians and obtained insightful results 
by studying, for example, systems of up to N e — 24 elec- 
trons at v = 1/3 ETT4261 . The infinite cylinder geometry has 
several additional numerical advantages, such as translation 
invariance both along and around the cylinder, the absence of 
curvature effects 11261 . and the absence of gapless edge ex- 
citations which slow down the convergence. In contrast to 
the standard bipartition of the torus lfT4l |27] 28 1, cutting the 
infinite cylinder gives the entanglement spectrum of a single 
edge. 

We show how to obtain several quantities that characterize 
topological order from the iDMRG simulations. The most im- 
portant step is to use iDMRG to systematically obtain the full 
set of degenerate ground states on an infinite cylinder, start- 
ing from a Hamiltonian. Their entanglement spectra encode 
the quantum dimensions d a lfi6l [TTl [29] and charges Q a of 
the quasiparticles. Topological quasiparticles on a cylinder 
appear as domain walls between pairs of degenerate ground- 
states, which we numerically optimize in order to obtain the 
energy and MPS of a localized quasiparticle of a chosen topo- 
logical charge. We apply the infinite cylinder technique to 
spin-polarized electrons at filling v = 1/3, 2/5 and 1/2 for a 
range of different Haldane pseudopotentials. 

Several topological properties of the FQH state manifest 
themselves on a torus rather than on an infinite cylinder. To 
construct the MPS for the ground-states of a torus with a 
'twist' in it, we cut out a segment of the infinite cylinder 
MPS and insert a "twist operator" before gluing the ends 
back together to form a ring, which is justified so long as 
the correlation length is small compared to the new circum- 
ference. By adiabatically changing the twist we can compute 
the non-Abelian Berry phase Ut associated with a modular 
T-transformation, which determines the topological spins 9 a 
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FIG. 1. Single particle states (p(x, y) in the lowest Landau level on 
the cylinder. L is the circumference of the cylinder and Ib is the 
magnetic length. 



of the quasiparticles, the central charge c of the edge theory, 
and the Hall viscosity r\ H of the bulk [30 31 1. Remarkably 
the result again depends only on the entanglement spectrum 
of the infinite cylinder. 

Model and method. We consider the cylinder geometry 
ll32T[34 | with a coordinate x running around the circumference 
of length L, and y running along the infinite length of the 
cylinder (see Fig. lib. The Landau gauge A = £g 2 (—y,0) 
conserves the x-momentum around the cylinder. The orbitals 
in the first Landau level take the form 
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where j g Z, with a density of one orbital per flux-quanta. 
Because each orbital is localized at yj — kjt 2 B , we treat the 
system as a ID chain. The most general two-body interaction 
allowed by translational symmetry can be specified by coeffi- 
cients V km , 



and MPS have been proven to be extremely successful in the 
simulation of gapped, one-dimensional systems because their 
ground-states can be expressed to very high accuracy by keep- 
ing a relatively small x even f° r infinite system size ll36l[39l . 
The iDMRG algorithm proceeds by iteratively minimizing 
E = {^\H\fy) within the space of iMPS. In the Supplemen- 
tary material we address several technical issues particular to 
QH iDMRG ED). 

After each variational optimization, the wavefunction is 
projected back to the original MPS form by keeping only the 
most important contributions in the Schmidt decomposition of 
each bond n, 
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The Schmidt states \a) L / R form orthonormal bases for the 
sites to the left and right of the bond n, and the Schmidt 
values are \ n[a . The entanglement entropy for this bipar- 
tition is directly obtained from the Schmidt values S n = 

We incorporate both particle number and momentum con- 
servation of | in its MPS representation and DMRG algo- 
rithm [41 1, which assume an important role in the MPS' sub- 
sequent analysis. The quantum numbers are defined to be 

C = C n = '^^(N n — v) (particle number), (6a) 

n n 

K = K n = n(N n — v) (momentum), (6b) 



H — ^3 ^ kmC n+m C n+k C n+m+k C n- 
n fe>|m| 



(2) where N n is the number operator at site n. The Schmidt states 
\ a ) lib, °f W a l so ^ ave a we ^ defined particle number C fi;Q , 



If the interactions are local in real space, they will be Gaus- 
sian localized with a range on the order of L/£b in the chain 
representation. For example, for the hard core Haldane pseu- 
dopotential V\ [32|, Vk m takes the form 



K 3 (k 2 



m 2 )e-^( k+m \ 



(3) 



To represent H at some fixed accuracy we must keep on the 
order of (L/£b) 2 terms in the Hamiltonian, and organizing 
these terms is one of several difficulties in the ID approach 
we address. 

For our numerical simulations, we use the iDMRG algo- 
rithm [35 1 which is based on the infinite MPS representation: 
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where fiWi™ are x x X matrices and \j n ), j n <= {0, 1} rep- 
resent the occupancy at orbital n. Assuming the state |^) is 
translationally invariant with a unit cell of length M, then we 
need only store M different tensors, i.e., = B^ l+M \ 
At filling v = p/q, M must be a multiple of q. DMRG 
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representing the total charge to the left of bond n. We will 
view On as a diagonal matrix acting in the set of Schmidt 
states (as for K). In what follows it is useful to define a 'bond 
expectation value', 



(C f , 



C n 



(7) 



which captures the expected value of the charge to the left 
of bond h (as for K). The Schmidt values and their quan- 
tum numbers, (A, C, K) n , constitute the 'orbital entangle- 
ment spectrum' (OES) of the bond. 

Topological order and quasiparticles. To interpret our re- 
sults we recall some basic facts about topological order on a 
cylinder ll42l . A topologically ordered state with m quasipar- 
ticle types (labeled by 'a') has an m-fold ground-state degen- 
eracy on an infinite cylinder. In addition, the chiral confor- 
mal field theory (CFT) describing the edge contains m scaling 
operators {</> Q } which insert a corresponding quasiparticle a 
near the edge. The Hilbert space of the edge H can be decom- 
posed into a direct sum of the Hilbert spaces associated with 
each quasiparticle, H = o /H a . Each space W a contains 
the spectrum of edge states with topological charge a. In ad- 
dition, the states of the low-lying entanglement spectrum are 
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FIG. 2. (Color online) (a) Entanglement entropy S of the v = 2/5 
state as a function of the circumference L, for a bipartition into 
two half-infinite cylinders. The inset shows the estimated TEE 
7 = Lj^ — S(L) converging to 7 = h log 5 (green line) for large L. 
(b) Estimate of the chiral central charge c/24, the topological spin h 
of the e/5 quasiparticle, and the Hall viscosity as expressed through 

the 'shift,' S = -j-rVn 1311 - Dashed lines show expected hierarchy 
values of 2/24, 1/5, 1/4 respectively. That h is identically correct is 
peculiar to abelian quasiparticles. 
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FIG. 3. (a) The relative energy difference between states in the 0110 
and 0101 sectors at v = 1/2 filling and L — 16£b- Using only 
the pseudopotentials Vi and V3 we adjust V3/V1 to tune between 
the CFL (small V3) and the MR state (large V3). Inset shows the 
real-space entanglement spectrum of the 0101 state plotted versus 
K. (b) Fixing a point V3/V1 = 0.4 in the MR phase, we increase L 
and measure the difference in entanglement entropies S0101 — 5*0110. 
The result is consistent with a quantum dimension d a — v2, the 
quantum dimension of the e/4 quasiparticle. 



in one-to-one correspondence with the states of the edge CFT 
ifTSl [TBI . With an entanglement cut running around the cir- 
cumference of the cylinder, we can choose a particular basis 
{|S a )} for the m-dimensional vector space of ground-states 
such that the low-lying entanglement spectrum of |H ) con- 
tains states only in W a . As each |H„) contains fewer states in 
the low-energy part of the entanglement spectrum than the full 
CFT, they have lower entanglement entropy and are referred 
to as the 'minimal entanglement states' (MESs) fl31 . Trans- 
lations permute the MESs, so the entanglement spectrum of a 
state depends on the bond n between adjacent orbitals where 
the cylinder is cut. Consequently each bond n is labeled by 
the sector a found in its entanglement spectrum; we therefore 
denote the corresponding bond by a [4-3 1 . 

We now outline the procedure for finding the full set of de- 
generate ground-states using iDMRG. On an infinite cylinder 
geometry the MESs are energy eigenstates, possibly with a 
degeneracy split exponentially in L fl4ll . Furthermore, for 
QH problems in a topological phase, initializing the iDMRG 
using a tensor product state of orbital configuration 'A' (for 
example, A = 010 gives | - ■ • 010010 • • ■)) will generate the 
set of distinct MESs states after iDMRG optimization. For 
many quantum Hall states, this is a consequence of MES' dis- 
tinct quantum numbers K, though we suspect the result is 
more general l40l . If we find optimized energies such that 
E\ 1 < E\ 2 , then the state derived from initial state A2 is 
rejected. The iDMRG is the numerical check that a given A 
leads to one of the m ground-states. The orbital string A asso- 
ciated to the MES |3 ) is called the state's 'root configuration' 
or 'pattern of zeros' B31I461 . 

The entanglement entropy of each MES |S ) scales with 
the circumference as S rj sL—j a , where j a = log \/^2 b d% — 
log (d a ) are the topological entanglement entropies (TEEs) 
and d a are the quantum dimensions of the quasiparticles 
|[T6l [T71 |29l . From the TEEs 7 a we can determine if we have 
the complete set of MESs [44|. 



As an example of an Abelian model (d a — 1 for all a), 
we consider v = 2/5 filling for Haldane pseudopotentials 
V3/V1 = 0.05, for which the seed |01010) and its translates 
were numerically determined to provide the 5 MESs. This has 
been considered before in an ED study with Coulomb interac- 
tions on a torus for circumferences up to L = 18£b H4], from 
which it was difficult to determine 7 accurately. We are able 
to go up to L = 23.5£b and obtain the results shown in Fig. [2] 
The estimated TEE is 7 » 0.83 (close to \ log 5 0.8047). 

As an example which has a non-Abelian phase, the filling 
v = 1/2 contains both the gapless composite Fermi liquid 
(CFL) phase and the gapped, non-Abelian MR phase Il42ll47l . 
Using the Haldane pseudopotentials V\ and V3, we tune be- 
tween the CFL (small V3/V1) and MR (intermediate V3/V1) 
phases. We start the DMRG either with the |01 10) configu- 
ration (which provides four states via translation) or |0101) 
(which provides two). We observe that in the suspected CFL 
phase the energies of the two sectors are split, while the MR 



phase is nearly 6-fold degenerate, as illustrated in Fig. 3(a) 



Fixing a point in the MR phase, the difference in the entangle- 
ment entropies of the MESs is S0101 — Sbno ~ 0.36, which 
implies that the Soioi state i s associated with a non-Abelian 
particle of quantum dimension 1.43 (close t o that of the e/4 
excitation: v2 s» 1.41), as illustrated in Fig. 3(b) This veri- 
fies the non-Abelian nature of the state. 

In the cylinder geometry, the topological quasiparticles ap- 
pear as domain walls between the degenerate ground-states. 
After numerically obtaining the infinite cylinder MESs, to 
generate a quasiparticle of type a in the vicinity of y = we 
use the matrices of the identity MPS |Si) for sites at y -C 
and the matrices of |H ) for y 3> 0. In the vicinity of y = 0, 
we insert a finite number of B-matrices which we numerically 
optimize in order to glue together the boundary conditions. 
For efficiency we work with quasiparticles of fixed momen- 
tum K; the resulting quasiparticles for the v — 1/3 state are 
illustrated in Fig. [4] To calculate the charge of the particle we 
measure Q a — lim e _j.o e~ e ' n 'C n , where C n is the on-site 
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FIG. 4. (a) Schematic of the MPS used to represent a quasiparticle 
'a', (b) The charge density p(y) of the numerically optimized ±e/3 
quasiparticles of the v = 1/3 state. The cylinder has circumference 
L = 16£b and the pseudopotentials up through Vi are chosen to 
approximate a dipolar r~ 3 interaction. 
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FIG. 5. (a) A torus with dimensions L x x L y . Taking the modular 
parameter r — s> r + 1 corresponds to a Dehn twist of the torus and 
generates the modular T-transformation on the set of ground states, 
(b) The wavefunction on the torus represented by a periodic MPS. 
The sites are labelled by n g Z + $5- due the flux threading in the 
^-direction. 



charge operator defined in Eq. ([§}. The resulting charge is in 
fact independent of the S-matrices used to glue together the 
MPS; Q a can be calculated exactly knowing only the OES of 
the infinite MPS |S„): 



a 2m(((C))-Ca) 



(8) 



where ((C)) is the average of (C fi ) over all bonds n in the 
MPS of |S a ) 1431 . This expression has been confirmed to be 
exact and reduces to the 'pattern of zeros' result in the limit of 
a thin cylinder [46|. 

Ground states on a torus and topological spin. Several 
quantities of interest, such as the phase's response to flux in- 
sertion and modular transformations, are defined on the torus 
geometry, so it is useful to convert the iMPS to a torus MPS. 
The torus is made by taking a cylinder of circumference L x 
and identifying points (x,y) to (x + T x L x ,y + L y ), where 
t = t x + iL y /L x is the modular parameter. We also allow 
fluxes Q x , $ y to thread through the two cycles. For a suitable 
definition of the iV$ = L x L y /2nt 2 B orbitals of the torus ll40l . 
we can construct wavefunctions on the torus from those on 
the infinite cylinder by taking a finite segment of the cylinder 
MPS and connecting the two edge auxiliary bonds together to 
form a ring [44], as illustrated in Fig.[5]l3). There is no need 
to reoptimize the B-matrices near the seam: locally the torus 
Hamiltonian is identical to that of a cylinder, and so long as 
L y is much greater than both the correlation length and the 
extent of the interactions, locally the periodic MPS will have 
correlations identical to the iMPS. The fluxes and modular pa- 
rameter can be accounted for by inserting a diagonal matrix G 
when connecting the two edge auxiliary bonds |40|, 
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where C and K are the conserved particle number and mo- 
mentum of the Schmidt states and N e = N<$,v is the total par- 
ticle number. The first factor enforces fermion statistics of the 
orbitals, the second factor adds a 2ttt x twist while connecting 



the ends of a cylinder to form the torus, and the final fac- 
tor arises from the flux Q y threaded through the y-direction. 
Via this construction we obtain the set of minimally-entangled 
states on a torus |H ) for arbitrary r and & x / y . By adiabati- 
cally varying these parameters, we obtain the associated Berry 
phases characteristic of the topological order in the system. 

The Berry phase Ut calculated as t. x goes from to 1 
||48l corresponds to a 'Dehn twist' of the torus as shown in 
Fig. [5ja). Ut is diagonal in the MES basis, and we expect 
it to contain two contributions, Eq. ( |10a| >. First, when acting 
on |H ), T causes an anyonic flux a to wind once around the 
x-cycle of the torus, generating a phase h a — St, where h a is 
the spin of a and c is the chiral central charge of the edge the- 
ory. Second, shearing the bulk introduces a phase due to the 
universal 'Hall viscosity' r\ H of the fluid 1301I3D . Ut can be 
calculated exactly by making use of the torus MPS [Eq. 
and we find that the result depends only on the infinite cylin- 
der OES, Eq. ( |10b| ). Equating the expected and exact results, 
for fermions 
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In Fig. [2] we use the v = 2/5 OES obtained from iDMRG 
to extract h, c and T) H , and find good agreement with the ex- 
pected values for the Abelian hierarchy state [49 1. For the MR 
phase, we obtain excellent results for the model Hamiltonian 
ll40ll . but when using only V\, V3 the measurement is highly 
sensitive to tunneling between the Pfaffian and anti-Pfaffian 
states present at the sizes studied ll50ll5TI . 

Conclusions. In this paper we showed how to numerically 
calculate several quantities which characterize topological or- 
der starting from a microscopic FQH Hamiltonian. The ap- 
proach consists of two key steps: (i) We find the MPS rep- 
resentation of the ground state using a variant of the DMRG 
algorithm, (ii) We derived expressions which allow us to ex- 
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tract the topological entanglement, quasiparticle charge, quan- 
tum dimensions, topological spins, quasiparticle charges, chi- 
ral central charge, and Hall viscosity of the phase from the 
MPS representation. As concrete examples, we applied our 
technique to spin-polarized electrons at filling v — 1/3, 2/5 
and 1/2 for a range of different Haldane pseudopotentials. 
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Independently, Cincio and Vidal have developed a similar 
technique for using DMRG to probe quasiparticles [52|. Also 
Tu, Zhang and Qi reported a similar method for extracting the 
topological spin from entanglement [53]. 
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ORBITAL n AND BOND n NOTATIONS 

In the subsequent analysis, sites are labeled by n, and bonds by n. The bonds n are understood to take on a numerical value 
indicating the location of the orbital, in units of 27 ^ B . As will be discussed in Eq. ( [37] i, when a flux Q x threads through the 
cylinder, n € Z + ^ . It is also convenient to understand n to take on a numerical value which is defined to be the average of 
the site locations to the bond's left and right. For example, the bond between sites n = and n = 1 is labeled by ?! = |. Bond 
n + M then refers to a bond M sites to the right of h. 

SYMMETRY PROPERTIES OF QUANTUM HALL IMPS 

Because the symmetry properties of MPS are an important part of the details that follow, we present the symmetries of the 
LLL problem, a review of U(l) symmetry conservation for MPS, and explain in more detail the notation (Q), ((Q)) notation 
used in the main text. 

Conserved quantities of the infinite cylinder 

In addition to the conservation of number, it is essential to implement conservation of momentum (also called 'center of 
mass'), both because of the tremendous reduction in numerical effort and for the ability to label the entanglement spectrum by 
momentum. At filling factor v, we define the charges (C, K) to be 

C = C n , C n = N n — v (particle number), (11a) 

n 

K = K n , K n = n(N n — v) (momentum), (lib) 

n 

where N n is the on-site number operator. Unlike the charge C, the momentum K is peculiar as it behaves nontrivially under a 
translation by n sites, T n , 

f- n KT n =K + nC (12) 

which must be accounted for when conserving the two charges in the DMRG algorithm. If the state has a unit cell of M, the 
matrices of the MPS are periodic in M ; however, due to Eq. ( fT2| ), the charges assigned to the Schmidt states are not. Under 
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translation, the charges of the Schmidt states on honda n, n + M are related by 

(Cn+M, Kn +M ) = (C n ,K n + MC n ~ M ((C))) . 

The value ((C)) is a constant of the MPS defined in the next section [Eq. (|20}]. 



(13) 



U(l) charge conservation for iMPS 

Let Qt = Yin Qn t> e an Abelian charge given by a sum of all single site terms Q n . Making a cut on bond n, we can decompose 
Qt = En<n Qn + E n >n Qn = Ql + Q R- For afinite size state \tp) with a Schmidt decomposition \ip) = J2 a ^ a \ a ) L \ a ) R 
on bond n, the Schmidt states must have definite charge, Ql \a) L — Q n . a \a) L . In the MPS representation, charge conservation 
is expressed through the corresponding constraint BTII541 



B afi -° 



where r/l denote the bonds to the right/left of site n. Pictorially, exponentiating the constraint implies 



(14) 



-□- 

1 



-□- 

A 



-A- 



(15) 



We can define a diagonal operator Q fl acting on the auxiliary bonds of the MPS with diagonal entries Q r - va . It is convenient 
to define a 'bond' expectation values of a Q n by 

(«fi) = E4^ (16) 

a 

where An are the Schmidt values on the bond n. In the finite case, this gives the expected charge to the left of the cut, (Qn) = 
(Ql). 

In the case of an iMPS, the necessary and sufficient condition for the iMPS to have definite charge is again that there exist bond 
operators Q fl such that Eq. ( p~4] > is satisfied. However, Eq. ([14} is clearly invariant under a uniform shift of the bond charges, 
Qn Qn + c ' so we can 't obviously interpret Q n as the physical charges of the Schmidt states. This ambiguity was absent in 
the finite case because the left-most 'bond' at the boundary can canonically be assigned charge 0. To resolve this ambiguity, we 
can explicitly calculate the charge of a Schmidt state using a regulator, Q a L — lim e _j.o Tlm<n e<Lm Qm \ct) L . Using Eq. ([14} we 
can rewrite this charge using that auxiliary operators Q 7Tl . Through an abuse of notation, we write Q rJl \ot) L as an operation on 
the state, by which we mean we insert Q m into the corresponding bond of the iMPS. We find 

Q aL \a) L = lim £ e £m (Q m+1/2 - Q TO _ 1/2 ) \a) L [Eq. ([14}] 

m<n 

W) L (17) 

Because e — > 0, the second contribution is independent of an arbitrary number of sites near the boundary n. Assuming the 
state has finite correlation length, far from n the result becomes independent of the choice of Schmidt state a. Making use of 
translation invariance (assuming a unit cell of length M) and the limit f->0we find 

HS& e E e£ *3*K = if E (Qrn)\a) L - (18) 
m<n Q<m<M 

where (Q m ) is again the bond expectation value of the infinite MPS. Hence the physical charge of the Schmidt state is 

QaL^Qn-^-jj E (Qrn) = Qn; a -((Q)), (19) 
0<m<M 



Qri 



lim e 

e->0 



E 
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where define the (()) notation: 

0<m<A/ 

The second term ((Q)) is an average of (Q l7l ) across the M sites of the unit cell. Using the freedom to shift Q — >• Q + c, we 
can always cancel the second term ((Q)) so that Q has its 'naive' interpretation as the charge of the left Schmidt state. We will 
call this the 'canonical' choice of bond charges. The state has fractional charges if the required constant c is a fraction of the 
elementary charge. 

NUMERICAL METHODS 

In the following we explain in detail three algorithmic issues particular to FQH DMRG on an infinite cylinder: construction 
of the matrix product operator (MPO) for the Hamiltonian, U(l) charge conservation for the momentum K around the cylinder, 
and ergodicity issues of the DMRG update. The MPO formulation of iDMRG used here is explained in Refs.l35landl55l 

MPO representation of the FQH Hamiltonian 

An MPO is a direct generalization of an MPS to the space of operators 15514571 . 

H= ]T ■■■®W£} a2 ®WW a3 ®--- . (21) 

0<cti<D 

Each is a matrix of operators acting on site n. The dimension D of the matrices depends on the Hamiltonian, increasing as 
longer range components are added. For the systems studied here, D ~ 100-300. An arbitrary two-body interaction in the LLL 
takes the form 

A 

H = U m nC i+ 2m+n C i+m+n C i+m C i + ^" C - 3 ^ mn '^i+2rn+n' l Pi+m+n' l f } i+m^i (^2) 

i 0<m.n i 0<m,n 

To represent the MPO exactly in the limit of an infinite cylinder would require taking D oo, but as the terms in the Hamiltonian 
of Eq. ( p2| ) decay like Gaussians when insertions become far apart, it is reasonable to truncate the Hamiltonian at some fixed 
accuracy by keeping only the largest terms. For simplicity we will use a cutoff to, n < A ~ L. 

To illustrate how to construct the MPO for Eq. ( p2) we view the MPO as a 'finite state machine' for constructing Hamiltonians 
l58l . Assuming translation invariance, to each index a of the matrix W a p we associate a state 'a' in a finite state machine, 




FIG. 6. Illustration of the finite state machine associated with the MPO representation of terms U mn in Eq. \22) . The machine begins in the 
node on the far left. When the first operator xpi is placed the machine enters the square grid. Each row of the grid corresponds to a value of m; 
each column, a value of n. The machine proceeds vertically up the left-most column until placing ip i+m , then proceeds through row m until 
ipl +m+n is placed with weight U mn , at which point it skips to the rightmost column. After descending down the rightmost column the final 
insertion ip i+ 2 m +n is placed, at which point the term is complete and the machine remains in the terminating node to the far right. The set of 
all routes from the far left to far right nodes generates precisely the terms U mn - A copy (with ip and ijy swapped) also exist (by not shown) to 
generate the Hermitian conjugate. 
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illustrated by a node in a graph. Each non-zero entry W a p is a transition probability in the finite state machine, illustrated with 
an edge. At the rtth step, if the machine makes the transition /3 — > a then the operator W a p is placed at site n. The machine is 
non-deterministic; if there are two possible transition out of the state /?, then the paths are taken in superposition, which generates 
the sum over all terms in the Hamiltonian. 

Assuming bosons and focusing on the particular contribution highlighted in Eq. ( |22] >, at each step the MPO will place one of 
1, tp or ijv . The MPO has a set of nodes which can be organized into a square grid, essentially in correspondence with the terms 
U mn , as explained and illustrated in Fig. [6] The rectangular nature of the graph leads to the scaling D ~ A 2 of the MPO. For 
fermions, the Jordan- Wigner string can be accounted for by replacing 1 with the string operator (— 1) F where appropriate. As 
for an MPS, U(l) charge conservation for both number and momenta can be implemented by assigning the appropriate charges 
(C, K) to the D indices of the auxiliary bond. 

Momentum conservation on an infinite cyiinder 

Knowing the charges transform as Eq. ( p"3j ), the algorithm proceeds as follows. We store the M i3-matrices of sites n — 
0, . . . , M — 1, and keep track of the auxiliary quantum numbers on the bonds to their right, n = |, §, . . . , M — |. For updates 
within the unit cell, charge conservation is implemented as usual. However, for an update acting on sites and M — 1, we must 
'translate' the charge data associated with site-0 to site- M using Eq. ( fT3j ). After updating B, the new charge data is translated 
back to site 0. 

Ergodicity of the iDMRG afgorithm 

The final peculiarity of applying iDMRG to the QH effect concerns the 'ergodicity' of the 2-site update. The standard iDMRG 
algorithm optimizes two neighboring B-matrices per step in order to avoid getting stuck in local minima of the energy landscape 
[20 , 35 . 55 1 . This has the added advantage that, unlike a 1 -site update, the bond dimension \ can grow during the simulation. For 
most Hamiltonians the 2-site update is sufficient to find the optimal state, even if the initial state is taken to be a product state. 
However, due to the additional constraint of momentum conservation for QH, starting from some particular state it is impossible 
for the 2-site update to generate amplitude in all possible configurations. The most naive explanation is that the smallest move 
available to the DMRG is a 'squeeze' involving 4 sites, though this picture is not quite exact. 

For the case of fermions, we have formalized and proven the following bound. Recall that because (C, K) are good quantum 
numbers, each Schmidt state a on bond fi can be assigned a definite quantum number (C n ; a , K n . a ). We define a combination 
P of these charges by 

— 1 -2 

Pn- a = K n . a - — C n . a - nCn;a- (assume ((C)) = 0) (23) 

P has been defined so as to be invariant under translation, unlike K. For the Laughlin states, which have an entanglement 
spectrum in one-to-one correspondence with a chiral CFT, P is precisely the total momentum of the CFT's oscillator modes 
lfl"8l . Let {P} be the set of P present on all bonds in the MPS, and let P m im Pmax be the minimum and maximum values they 
take before beginning of DMRG. We have proved that using a 2-site DMRG update {P} always remains bounded by P m i n and 
Pmax- Hence the entanglement spectrum will appear to have a momentum cutoff set by the initial state, and the 2-site update 
will fail to find a variationally optimal state. For example, if we use the exact Laughlin state as the seed for DMRG (which 
has P m in = 0, as the model state is purely 'squeezed'), the 2-site update will not arrive at the ground state of the Coulomb 
Hamiltonian, which has P m i n < 0. 

Though not stated in these terms, to our knowledge this ergodicity problem was previously dealt with via two methods. One 
approach initialized the DMRG with a large spectrum of random initial fluctuations, and ensured by hand that the DMRG update 
preserves several states in each charge sector (e.g. Ref. 11251 ). The algorithm was nevertheless observed to get stuck for several 
sweeps at a time, but did converge. In this approach care is required to ensure the initial state supplies adequate fluctuations 
Pmin, Pmax> an d there may be additional more subtle restrictions missed by this bound. A second approach uses a two site 
update with 'density matrix corrections' (e.g. Ref. l23l ). 

The solution we adopt here is to generalize the DMRG to a ri-site update, optimizing n sites as a time. In the MPS/MPO 
formulation of DMRG we implement a 2n-site update by grouping n sites together into a single site of dimension d = 2™ and 
then perform the usual 2-site update on these grouped sites. For n = 2, for example, we simply contract 2 adjacent S-matrices 
of the MPS, taking B^B^ 1 -> B^ jojl , and likewise for the V^-matrices of the MPO. As the complexity of the DMRG 
update scales as d 3 , this does come at a cost, though it is partially offset by the increased speed of convergence. We have not seen 
the algorithm get stuck while using a sufficiently expanded update. Though we have not proved it, we believe that ag + 3-site 
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update will be ergodic for filling fractions ^, 1 — |. For more complicated fractions, such as |, we have checked ergodicity by 
trial and error; for instance, a 10-site update arrived at the same final state as a 6-site update, so the latter was used in the reported 
simulations. 

The main advantage of previous approaches is the decreased memory required, which quickly becomes a limitation for the 
6-site update when \ ^ 4000. However, generating a random iMPS with sufficient fluctuations, conserves K, and resides in the 
desired topological sector is quite subtle. An intermediate approach is to use a 2n-site update for the initial steps to generate the 
necessary fluctuations, then switch to a 2-site update until convergence is reached. 

Obtaining the full set of minimal entangled ground states {|H a }} from iDMRG 

There are two issues when constructing the MES basis, a) Does the DMRG converge to MESs, or is it favorable for it to 
converge to superpositions of them? b) If it does produce MES states, how do you initialize the DMRG in order to obtain all of 
them? 

As to a), it has been shown that the MES basis is in fact the eigenstate basis on an infinite cylinder, with energy densities that 
have an exponentially small splitting at finite circumference L [44]. As there is no energetic reason for the DMRG to produce 
superpositions of the MES (and at finite x, there is in fact a finite entanglement bias to produce MES), we expect the DMRG to 
produce MES. 

As to b), we first consider the role of the momentum K per unit cell. The MESs are eigenstates of the momentum K: if two 
infinite cylinder ground states with different momenta per unit cell are added in superposition, then the entanglement entropy 
must increase as their Schmidt states are respectively orthogonal. The DMRG preserves K; so if the DMRG is initialized using 
a state in a momentum sector K that contains an MES, then the DMRG will produce an MES; if the sector does not contain an 
MES, the optimized energy will observed to be higher than those of sectors that do. The first step, then, is choose an orbital 
configuration A of the desired K (such as A = 010), initialize the iDMRG with the corresponding \ = 1 MPS |A) by updating 
the iDMRG 'environments' without further optimizing the state, and then run iDMRG to obtain an optimized state and energy 
E\. Repeating for orbital configuration A of different K, we compare the energies E\ to determine which K sectors contain a 
MES, rejecting those sectors such that E\ is not minimal within some tolerance set by the exponentially small splitting due to 
L. 

For many of the expected phases (such at the Laughlin and hierarchy states), the ground-states are uniquely distinguished by 
K, so all MES will be obtained by the procedure just outlined. For certain cases, however, several of the MES have the same K, 
for instance those corresponding to root configurations OHIO and 10101 of the k — 3 Read-Rezayi (RR) phase at v = I l59l . 
If we initialize the DMRG with 10101, we might worry that it will tunnel into the 01 1 10 MES during the iDMRG, either due to 
the exponentially small splitting of the physical energies, or because the latter state has d a > 1 and hence higher entanglement. 

We argue that the iDMRG, if initialized with an approximation of one MES 6, will not tunnel into a different MES a. Consider 
the following three energy scales: 1) the exponentially small, but physical, splitting between the ground-state energy per site, 
Eo-,ab = Eo- a — Bo;6'> 2) the difference in the DMRG truncation error E x:a b per site (inherent to the MPS representation), each 
of which can be made arbitrarily small for sufficient x, but may still differ if d a ^ d^; 3) the gap A c for inserting a quasiparticle 
of the type 'c' that would arise at a domain wall between the two states a, b. 

If ^0;a6 + E x - a i, < 0, heuristically the iDMRG may prefer to find the a state. In this scenario, if we initialize the iDMRG with 
the b state, it can 'tunnel' into the a state by inserting a c, c pair within a correlation length £ of the sites being updated. As the 
iDMRG proceeds, the state 'grows' by repeatedly inserting new sites at the center of the chain. The c, c pair then get successively 
pushed out to the left/right of the chain, leaving the a type GS in the central region, which the state eventually converges to. 
Effectively, a c/c pair has been drawn out to the edges of the 'infinite' cylinder, thus tunneling between ground-states, a problem 
the geometry is supposed to avoid. 

However, we do not expect this to happen for simple energetic reasons. The cost to tunnel into the a state at a given update 
step is 

E = A c - (£E . ab + E x . ab ) (24) 

where £ is the correlation length of the system (in sites). Eo- a i, is exponentially small at large L, and E x should be very small if 
sufficient \ is used, while the quasiparticle gap A c remains finite. Hence the energy of the quasiparticle provides an energetic 
barrier for the iDMRG to tunnel between the MES. 

If the root configurations A are close enough to the desired MES for the purposes of the above argument, then by initializing 
the DMRG with different A of the same K, the DMRG should produce the corresponding orthogonal MES. Testing successively 
more complicated orbitals, we can check if we have obtained a full set by summing the quantum dimensions of the states 
accepted thus far P4l . We have not verified if this proves to be the case for a non-trivial case such at the k = 3 RR state. 
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If the root configurations A are not sufficiently close to the MES for the purposes of the above argument, it is also possible to 
run iDMRG while including a bias against the MES obtained so far in order to find the additional MES. 



QUASIPARTICLE CHARGES 

In the main text we claimed that the quasiparticle charges Q a are determined entirely by the entanglement spectrum of the 
iMPS |S ), which we demonstrate here in detail. As discussed, we suppose the MPS takes the form of |3j.) for y < and |3 n ) 
for y > 0. The most general form such a state can take is 



E 

a/3 



a) L ®\a(i) c \P) R (25) 



where \a) L are the left Schmidt states of |3i), \/3) R are the right Schmidt states of |3 a ), and \a/3) c is an arbitrary set of states 
in the central 'gluing' region. Without loss of generality, we suppose that the the gluing region has a length which is a multiple 
of q, n = {0, 1, • • • , ql — 1} for I e Z. The boundary bonds til = 1, n>R = a are indexed by a/f3 respectively. Schematically, 
we should think of the low-lying Schmidt states on the left/right bond as being in the CFT sector Vi /V a respectively. 

We exploit three basic facts. First, because the central region contains a multiple of q sites, the charge C of \aj3) c must be 
a multiple of the electron charge (which has been chosen to be 1) for all a, (3. Second, the charges of the left Schmidt states 
are (Cj — ((C)) j), which differ from each other only by multiples of the electron charge. Likewise, the charges of the right 
Schmidt states are — (C a — ((C)) ), which differ from each other only by multiples of the electron charge. Hence the charge of 
the quasiparticle, Q a , satisfies 

e 2mQ a = e 27ti[(C I -{(C)) 1 )-(C a -((C))J.] (26) 

Again, we have taken advantage of the fact that C rl is a constant modulo 1, the charge of an electron. Finally, because charge 
conjugation C acts as spatial inversion on the LLL orbitals, we must have e 2 ^ l ^ Cl ~^ c ^ ^ = 1, As claimed, 

e 2irtQ a = e 2™«(C)) -C a )_ (27) 

When cutting a MES on any bond, C mod 1 is single-valued, so the equation has no ambiguity. 

Note that we have implicitly assumed the sector 1 appeared in the OES - but this need not be the case, as for instance in the 
v = h bosonic Laughlin state. If we naively apply the above formula to v = ~ Laughlin, we find particles of charge ±e/4, 



rather than 0, e/2. This issues is clarified using an alternate derivation of the charge via flux matrices 



MATRIX PRODUCT STATES ON A TORUS 



Here we explain in detail the procedure to take an infinite cylinder iMPS to a torus MPS. Our approach is closely related to 
that of Ref.04] but with the additional complication of having twisted boundary conditions. 



From an infinite chain to a periodic one 

We first step back and show how to convert the MPS of any gapped, infinite chain to a periodic chain, both for the trivial case 
when there is no 'twist,' and then in the presence of a twist generated by a symmetry. The construction in the first case, which 
we explain for completeness, is obvious: we cut out a segment of the iMPS and reconnect the two dangling bonds to form a 
ring. For simplicity in what follows we will assume bosonic chains, and introduce the correction for Fermionic chains due to the 
Jordan- Wigner string at a later point. 

Consider two systems, an infinite chain with a unit cell of N, and a periodic chain of length N, The sites of the chains are 
labeled by 'n' (with n ~ n+N in the periodic case). Restricting to local Hamiltonians, we can decompose H — ^2 n W n ' , where 
each is localized around site n over a length <C N (for bosons, for example, these are the terms ff W = b\ l+1 b n + h.c, 
though could extend over many sites). In the infinite case, translation symmetry implies T N t~ n = H^ n+N \ T being the 
translation operator. 

The energetics of the state are determined by the reduced density matrices of the system, E = ^ Tr(pM^H), The pt n l 
are reduced density matrices in a region around n large enough to include all the sites affected by W n \ If the ffW of the 
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finite and infinite chains are identical, then to find the ground state of the periodic system it will be sufficient to reproduce the 
local density matrices of the infinite system. To do so, we first cut out a segment of the iMPS with two dangling bonds, 
* a/3 = J3[°] J3W ■ • • Bl 1 *- 1 ], or in the pictorial representation of MPS: 



a- 



-□- 



-□- 



-□- 



-□- 



(28) 



1 N-2 N-l 

We then connect (trace over) the dangling bonds to form a ring MPS: 



Tr[*] 



□ - 



-□- 



-□- 



-□ 



N-2 



N - 1 







(29) 



There is no need to further optimize the MPS at the 'seam,' because when the complement to the region of is large compared 
to the correlation length, N — 3> £, then up to corrections of order e~ N ^ the of the periodic MPS are identical to those 
of the iMPS. 



Twists 



We now introduce a twist generated by a local symmetry, Q = J2 n Q^ n \ assuming each is individually symmetric under 
Q and that t n Q^T~ n = Q^ n+N ^. We use a unitary 'twist' operator G to introduce a twist every N sites of infinite chain, 
which generates a twisted Hamiltonian: 



N-l A[n + bN] 



H oo {6) = G{0)H oo G{-e)=Y,H [n] {6). 



(30) 



The chain remains translation invariant, f N H^(9)f- N = H^ n+N ^{9), because Q is a symmetry. 

We define the twisted Hamiltonian of the periodic chain to be H a {9) — J2 n =o (^) w ^ W ta ken fi" om the infinite 
chain. In the bosonic case, for example, there is a single link with a twist, b\b-\e %e + h.c. Note that H o (0) is not unitarily 
related to H o (0) except when 6 is a multiple of a 'flux quantum' $, which we can assume is $ = 2ir . What is the ground state 
\6) o ofH o (0)7 

By construction, the Hamiltonians H ,Hoo remain locally identical, so given the iMPS for |0) , we can utilize the gluing 
trick of Fig. d29| already justified. As |0) = G(Q) |0) , with 10)^ the untwisted ground state, the desired iMPS is 



-□- 



-□- 



-□- 



-□- 



-□- 



l^)oo= 

-1 1 2 

and so on throughout the chain. Using the conservation rule of \\5) , we can rewrite the above as 

-%9Q 

... □ u— 

l*>c 



-□- 



-□- 



-□- 



-1 







To obtain \0) o , we again cut out a segment and glue, 



□- 



G 

-□ • □- 



N-2 N-l 



-□ 



(31) 



(32) 



(33) 



where G = e 
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Fermions 



One additional modification must be made for fermions, due to the Jordan Wigner string. MPS for fermionic chains are 
always expressed through the occupation of bosonic operators er+ = (— l)5-^<» Ni ^ n , where iVj is the occupation at site i. In 
other words, the B matrix at site n generates the state according to B^™ (cr^)- 7 " |0). On a periodic chain, the fermionic operators 
ipn ~ ipn+N are periodic, but because of the string the er+ are not. If N is the bond operator corresponding to number, we find 
the string can be accounted for via 



G = v ( NF -V e e- i8 Q, r? = ±forB/F 



(34) 



where N F is the total particle number and r\ — ±1 for bosons and fermions respectively. 

In the following sections, we prese nt the det ails for this construction pertaining to quantum Hall systems. First we de scribe 
the single-particle basis on an infinite cylinder by solving the Hamiltonian H Q = |(p + A) 2 , then adapt the basis to a torus 
We then show that in the orbital basis, the cylinder and torus Hamiltonians H^, H Q , satisfy the criteria just discussed, so we can 



obtain the torus MPS from the cylinder MPS. 



The cylinder geometry 

Consider a cylinder finite in the x-direction with circumference L x , infinite in the y-direction, with the following boundary 
condition and Hamiltonian: 

i>{x,y) = iP(x-L x ,y)e i **, (35) 
tf = i(p + A) 2 , (36) 

where $ x is the flux threading the x-cycle and p = — iV. In the Landau gauge A = £g(— y,0), the x-momentum k is a 
conserved quantity. In that eigenbasis, the wave functions ifin l (x, y) — (x, y|yVi yl ) in the lowest Landau levels (LLL) are 

^\x,y) = — e ^ e - (V i'l , withfc n = n^, y n = k n £%, neZ+f^. (37) 

We see that the centers of the wave functions y n depends on $ x ; pumping a flux through x-loop will shift all the orbitals by unit 
distance Ay = 2irl 2 B /L x . Although the phases of these orbitals are arbitrary, it is convenient to resolve the ambiguity via the 
translation operator: 



T y — exp 



(38) 



which shifts the wavefunctions by Ay. We have chosen the phases of the orbitals such that |<^j y +1 ) = T y \(p„ ). 

Note. We emphasize that the flux <& x is accounted for by labeling sites as n € Z + This carries over to the definition of 
the momentum: K n = nC n includes the fractional part 5*. This definition of K n is an important detail for the formulas that 
follow. 



The torus geometry 

To go from a cylinder to a torus, we identify points (x, y) with (x — r x L x ,y — L y ) as follows 

i/j(x, y) = i>{x - t x L x , y - L y )e U B 2 L y ^ y (39) 

or equivalently ?p(x, y) = ip(x + T x L x ,y + L y )e~ ie B ^(i+t.l,)-!*,^ The phase i€g 2 L a a; is necessary since the gauge potential 
A is not periodic when taking y — > y + L y , the difference being A(x — r x L x ,y — L y ) — A(x, y) = V(£g 2 L y x). Note that for 
the above to be well-defined, we need an integer number of fluxes in the torus: 

^=»,€Z. (40, 
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The fluxes $ y ) parameterizes a set of Hamiltonians and corresponding Hilbert space of ground states. 

To figure out the eigenstates for the torus, we sum over combinations of </^ yl on the cylinder such that they respect the boundary 
condition above. The single-particle wave functions ip n (x, y) = (x, y\(p n ) f° r the torus are 

b 

We will use as our basis the orbitals < n < iV$. 

MPS for twisted tori 

Since the flux $ a corresponds to a twist generated by the charge C, while the modular parameter t x corresponds to a twist 
generated by the momentum operator K, it is natural to suppose the orbital MPS can be obtained using a twist generated by 
2itt x K — $> y C. We prove this intuition is correct, but the surprise is that the ground state remains exact (up to corrections 
Q- L y/£) even though the twist is performed in orbital space, rather than real space. 

We first determine the form of the orbital Hamiltonian of the torus, H Q (8), in the basis of Eq. ( pTTj i. We assume the Hamiltonian 
arises from products of the electron density operators p(x), is translation invariant in both x and y, and that any interactions, 
such as p(x)V(x — y)p(x), are short ranged in comparison to the length L y of the torus. As we will eventually take L y — > oo, 
this is not really a restriction. 

We make use of the interaction overlap integrals for the cylinder orbitals <p^}, which we denote Vs ni \. The two body term, 
for instance, is encoded in a term V nin2n3rii . Vi ni \ is non-zero only if ^ n, = 0, due to momentum conservation around the 
cylinder, and V/ ni _|_H = V{rii} due magnetic translation invariance along the cylinder. Furthermore, V} ni \ decays when the 
separation between indices are such that L x /£b <C |«i — rtj\. Since L x /£b <C N<&, to exponentially good accuracy we can 
assume Vs ni \ — if any two rii differ by Nq> sites. 

For notational simplicity we will illustrate the calculation for the 1-body term, which generalizes in an obvious fashion. The 
1-body Hamiltonian is 

H= £ ]T K 1+bl ^,n 2+b2 iv.e-^"- bl ) + ^"- fc2 ) C t i c„ 2 (42) 

0<ni<N<s, 0<6i<l 

f(n, b) = -2iTiN ^ b ~ l K x ~ 2mbnT x + ib^ y (43) 

The summand is invariant under taking all bi —> bi + 1, so we can safely restrict to terms such that bi € {0, 1} with at least one 
of the b, t — ( V vanishes if two indices are iV$ apart). This leads to two types of terms: the 'bulk' terms, in which all bi = 0, 
and the 'seam' terms, in which some bi = and some bi = 1. The seam terms arise when site at the beginning and end of the 
unit cell interact. 

Since f(n, 0) = 0, in the bulk the local Hamiltonians H^(6) = H^(0) are identical to the infinite cylinder Hamiltonians. 
Along the seam the (6) acquire phases f(rii, 1) = —2i:iT x ni + i$ y for each orbital with bi = 1. For example, if iV$ = 10, 
in the interaction Vc^nCgCi — > Vg^iiCgCie^ 1,1 -'. 

Consequently the torus should be identical to the infinite cylinder in the bulk, but differ by phases f(rii, 1) near the 
seam. To account for the phase f(rii, 1) = —2-KiT x rii + we use the same conserved quantity manipulations as we did in 
Eq. ( |32| ), and find the correct twist operator G is indeed 

G = rj^"- 1 ^ exp [-2niT x K + i%G] . (44) 

i] = ±1 accounts for bosons and fermions respectively, as discussed for Eq. ( |34) , There is one subtlety we must emphasize: 
recall that K is not periodic, due to Eq. ( fT2) , so the quantum numbers K a according to the right bond of the last site and the left 
bond of the first site differ by iV$C. In the above form we assume K uses the quantum numbers according to the bond before 
the first site, at n = — |, not after the last site. The first factor rf N reflects the fermion statistics among the orbitals. 

COMPUTATION OF FLUX AND MODULAR MATRICES VIA THE ENTANGLEMENT SPECTRUM 

The flux matrices T XfV gives the action of threading a 2-k flux in the x, y-loop, while the modular T-matrix gives the action of a 
Dehn twist. To derive their expression, we compute the (non-Abelian) Berry phase from the adiabatic changes $ x iV —> Q x , y + 2ir 



16 



and t x — >• t x + 1, respectively. In each case a parameter k is varied and the Berry phase U is a m x m matrix which is a product 
of two pieces: 

V = W foe* , (45) 

with A and W being matrices defined as 

A ab (n) = (3(,(/c) | -i^ | 3 a (/c)>, J¥ afc = (3 6 (k/) | S„ («,•)) . (46) 

As k is varied, |S a (/c)) denotes the set of ground states for the boundary conditions specified by k; A is the Berry connection 
and W is the overlap between the initial and final set of MES. While the result is independent of the choice of basis, in the MPS 
construction |3 («;j)) remains the MES basis. 

The action of T x / y ,T are characteristic of the topological phase and robust to perturbations. The results presented are 
summarized in Tab. [I] Although the flux matrices and topological spins are formalized for a state in the torus geometry, we will 
show that these quantities only depend on the OES of an infinite cylinder. In particular, knowing the entanglement spectrum of 
each MES |H ) allows us to determine the charge Q a and topological spin h a of the quasiparticle a. 

Note. There is a subtle point in the interpretation of the resulting equations. The charge Q a and the spin h a we will sub- 
sequently derive are those of the sector H a that would arise in the real space entanglement spectra if we were to cut the system 
at y — 0. The entanglement spectra appearing at y = depend on $ x . It is most natural to choose $> x such that a = 1 can 
appear on the cut at y = 0: the resulting set of 'a' are then those that would appear on the plane. We find that for fermions, this 
occurs for $> x = tt, while for bosons, this occurs for Q x = 0. Hence it is most natural to set <J> X = for bosons, and & x = w for 
fermions in the following equations. 



K parameter Berry phase matrix 


Physical observable 


Formula 


®x Tx 


translation structure 


Eq. {58} 




charge 


Eq. {56} 


Tx T topolo; 


jical spin, Hall viscosity, central charj 


;e Eq. {6T} 



TABLE I. Summary of the non-Abelian phases from pumping <fr x , <E> H and r. 



Computing the Berry phase 

We have already determined the many body state |H„(«)) = i &j a j 1 ...{ K ) \lo,ji ■ * • ; At), so it is in principle a mechanical matter 
to calculate the Berry connection. The orbital wave function ^ a is given by the periodic MPS with a twist [Eq. {33)], and j n is 
the occupation number of orbital ip n defined in Eq. {4T) . Our approach is as follows: we first derive the general structure of the 
Berry phase, by decomposing the result into the 'wavefunction' and 'orbital' parts. We then report the results for the constituent 
terms for T x i y , T. 

• A. The connection A a b(K) can be separated into two components, one from the changing twist in the MPS and one 
from the changing orbitals: 

The first term of Eq. {47) involves the Berry connection for 1 i ,a , which depends on At through the twist operator G. In the 
limit of L y 3> £, the required overlap reduces to a 'bond expectation value' as defined in Eq. {T6) , 

4?=-^ 6 (^f ). (48) 

The expectation value is taken on the bond where G sits, to the left site n = fe. The second term of Eq. {47) involves the 
Berry connection between orbitals: 

\(<p)r„-\— ;x . /at \ i« i,„ \ a(<p)i\T\= f KJ Aif), 



A%\K) = -i5 ab J2{N n ){Vn\d K \ V > n ), 6^{N n )= A^ b \K)dn (49) 

n jK i 

where (N n ) is the average occupation at site n. 
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• W. Likewise, the overlap W a b also contains a twist and orbital component, 

w ab = J2 %(«/) 

We distinguish between two cases: for and T, the orbitals return to themselves (), giving 

W <* = E %("/)*{*}(«<) = Sab^^- (for T y , T). 
For J^, the orbitals are translated by 1 site, and we find 

W ab = {E b \f- l \E a )^^ (for J*). 

l*\K=Kf 

Because of this decomposition, the total Berry phase takes the form 



U a b = W ab exp 



(50) 



(51) 



(52) 



(53) 



Below we summarize the results for each case, giving an explicit formula for F x , T y and T ■ 



Orbital Berry Phase 



A technical difficulty arises when we try to compute the orbital Berry phase 0« , due to the changing boundary conditions. 
To remove all ambiguity, we follow the approach of Ref. [48 1 by working in a coordinate system X, Y with fixed boundary 
conditions (which determines the Hilbert space) and vary the metric and gauge potential, which determines the Hamiltonian. In 
this approach, the boundary conditions, metric, and vector potential read 



il>(X,Y) =iP(X + L x ,Y), 

4>(X, Y) = fax, Y + L y ) e- a B 2 L v x 

( ^ Tel 



with t v = L y /L x , 



L, B ' L 



y L y J 



(54a) 
(54b) 

(54c) 
(54d) 



The odd looking term 7rT g JV '»' g Ay exists to counteract the N$/2 flux quanta inserted as t x increases by 1. 

The coordinate system (x, y) and boundary conditions of Eq. ([39]) are unitarily related to (X, Y) through a change of coordi- 
nates and gauge transformation. While we omit the details of the computation, all the orbital overlaps and Berry phases can be 
unambiguously defined by transforming the orbitals to the (X, Y) coordinates. 

We also note that under T, examining the single particle orbitals discussed above shows that the fluxes transform as 
(Q x , <f> y ) —> (<& x , $ a + Q x ), so Ut does not take the system back to itself unless $ x — 0. We return to this issue when 
discussing modular transformations. 



Flux matrices T x / y : quasiparticle charge and Hall conductance 

• J- y . The flux matrix T y describe how the MES |3 n ) transform as a flux quanta is threaded through the y-loop. Letting 
k = & y , the pieces in Eq. |53) are as follows, 



l v) = -2tt 



n 



= (C) 



W ah = 6, 



abe 



-2-rriCa 



(55) 
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The number operator N may be rewritten in terms of bond operators, e.g. N n = C n + v = C w+1 ; 2 — ^ n -i/2 Vk With 
a bit of algebraic manipulation, the Berry phase can be written as 

T y = exp [2m({(0)) - 0) + ii/($ a - it)] (56) 

We see that our formula for T y is similar to that for the charge e 27rl< 2» in Eq. ( p7| ). Note that is diagonal in the basis of 
the MES | S ). Physically, this is because each quasiparticle type has a well-defined charge, modulo the electron charge. 
In Eq. ( |56] > we explicitly include the Q x dependence, whereas in Eq. ( |27| i we implicitly assumed $> x = it, which is most 
natural for fermions. Since § x translates the state (see below), the term v$ x encodes the Su-Schrieffer counting. 

• T x . Inserting a flux quanta through the x-loop will in effect translate the MPS by one unit cell. We can see that the 
orbitals momenta, in units of fp, are quantized and of the form n = Z + fe-, and hence as Q x adiabatically increases by 

2ir, the orbitals ip n will evolve in to (p n +i- It is easy to see that T x must proportional to the translation operator T y . For 
k = Q x , the pieces of Eq. ([53]) are 

6^=2ttt x ~, A^ = -t x (C), W ab = 8 a . b+l e 2 ^- (57) 

(6+1 refers to the MES which is b translated by one site). The flux matrix takes the form 

T x = 5 a ^ b+1 Fy T - exp [iu(-$ y - n)] . (58) 

In general, the flux matrices must satisfy the algebra 

J~ X J~ y 6 yF X . (59) 



Dehn twist T: topological spin, central charge, and Hall viscosity 



Varying k = t x from to 1, we get the Berry phase Ut- Note that under T, the fluxes transform as T : ($ x , ®y) 
(Q x , $ a + Q x ). Thus we expect that at $ x = 0, Ut should be unambiguously defined, while at <f> x — ir, only U T should be 
unambiguously defined, a point we will return to later. 

For the / th Landau level, the individual parts of Eq. ( f53| ) are 



2N* 



4L„ 



=-(K), 



Wab = S ab e 



2itiK & 



(60) 



(The lowest Landau level is labeled by I = 0.) We remind the reader that for & x ^ 0, the site index n and K n include a fractional 



part see Eq. 



where 
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We note that the first term of 



U T -ab = S ab exp <^ 2m 



K-((H)) 



has been derived by Wen and Wang in Ref.06] Combining them, we have 

'$2 $, 1\ (2l + l)v 



2 V4tt 2 2tt 6 



16tt 2 



Li 



(61) 



H n = K n -nC n + n((C)). 



(62) 



In light of Eq. (f2"3]> (where ((C)) = 0), we can write H = P + ^-C and interpret H as the "energy operator" for auxiliary states. 



H is convenient as it is invariant under translation by q sites. As mentioned, at & x = 0, Ut is unambiguously defined because 
K mod 1 is single-valued on any bond (for a MES) and so Eq. ( |6Tj ) is well-defined. For <fr x — it, however, only K mod \ is 



single-valued which implies that only U T is well-defined (cf. section on modular transformations!). 

If we use a convention where ((C)) = 0, set & x = it (which is most natural for fermions) and restrict to the lowest Landau 
level, then Eq. ( |6"T] i simplifies to 



Ur-ab = exp ^2iri 



K - ((K - nC n )) - ^ - 

24 1&TT Z 



Li 



(When I = 0, ((C)) = 0, and <S> X = tt.) 



(63) 



This is the form shown in the main text. 




FIG. 7. Various quantities characterizing the topological Moore-Read phase at v = 1/2. The Berry phase Ut arising from a Dehn twist is 
acquired from the model Moore-Read wave function via Eq. (63} and h, c, rj H are extracted via Eq. l |65[ ) as function of circumference L — L x . 
(a) The argument of the phase Ut plot vs. L 2 /£%. For large L, the argument becomes linear in L . (b) The topological spin h of the a± 
quasiparticle, extracted from the ratio of Ut between the a± and 1 ground states, exp 2nih a± is given as the ratio of the topological spin 
8a ± /dt- (c) The Hall viscosity r] H extracted by fitting to the form Utm = 9i exp(— ^^L 2 ) for the identity sector. The data is presented as 
the 'shift' S = (2nl 2 B /i>){A'q H /K). (d) The chiral central charge extracted from Ut-,u, assuming S = 3 for the Moore-Read state, (a)-(d) In 
all cases L/Ib must be sufficiently large for the result to be reliably extracted from the entanglement spectrum. 



Discussion. In order to interpret Ut, we picture the MES as a torus with a topological flux a winding around the y-cycle. 
We cut the torus at y = 0, shear the segment of cylinder evenly throughout the bulk so that y = remains fixed while y = L y 
rotates by L x , then reglue the ends. Since an anyonic flux a terminates at the edge, there is a quasiparticle a on the edge y = L y 
which moves once around the circumference under the shear. The quasiparticle has momentum j^(h a — c/24) (if the edge is 
not chiral, this should be understood as h a — h a — (c — c)/24), hence traversing a distance L x generates a phase h a — c/24. 
The bulk itself shears as well, with the strain changing as kf-dr x . The finite 'Hall viscosity' r\ H results in a phase per unit area 
proportional to the changing strain, 



Ibulk 



(L x L y )?] H I ^dT x = h 1 rj H L 2 x . (64) 



=0 



Together, the expected result is 



U T ;ab = 5ab exp 



2m Ihn-^r- ^rl? 



(ha 



24 2nh 



(65) 



The Hall viscosity is known to be related to the 'shift' S on a sphere via r\ H = | 2 ^ gi § l3D . 

To extract the quantities independently, we can first fit the quadratic part to extract T) H and isolate the constant part h a — c/24. 
Comparing h a — c/24 across MES, we can identify the 1 sector because ht = minimizes this constant (strictly speaking, we 
only measure h a — c/24 modulo either 1 or 1/2, for & x — 0, ir.) Once 1 is determined, we obtain c, and the remaining h a 's can 
be read off directly. 

In Ref. 46 , the differences h a — hb were related to the 'pattern of zeros,' which can be understood as the L — > limit in 
which the state becomes a \ = 1 tensor product. In the L — > limit, where \ — 1> we have checked that the differences 
h a — hi, inferred from Eq. [6T|reproduces Wen and Wang's [46 1 result. However, in the limit L — > oo, the ratio h a — hb computed 
via Eq. ( |6T| matches that of Ref. 06| only when the quasiparticles a, b are related by attaching fractional fluxes (in the MES 
language, the states are related by translation). Hence in the (single-Landau level) Abelian case, where all MES are related by 
flux attachment, all h a — hb can be recovered in the limit L 0. But in the non-Abelian case, where not all MES are related by 
translation, only the L^oo limit gives the correct result, so the result of Ref. 46 does not appear to be correct. 

As an example, consider the Moore -Read state at v = |. The MR state contains two non-Abelian quasiparticle '<r + ' and 
'<7_', with charges ±| and both with quantum dimension d a = \2. Its topological spin is expected to be h a± — j^Q 2 + h a — 
+ = |. The edge of the MR phase consists of a free boson and Majorana mode with combined chiral central charge of 
1 + | = |. Ia Pig. [7J we plot the values of h, c and r\ H extracted by fitting Ut of the model MR wavefunctions B2l to Eq. ( |65j ) 
at various L x . 
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Modular transformations 



The modular transformations are the set of affine transformations from the torus to itself. For example, the 'T" transformation 
corresponds to a Dehn twist sending r — » r + 1 (the same as t x — > + 1). The '5' transformation rotates the torus sending 
t — > — 1/r. (When = 0, this corresponds to a 7r/2 rotation swapping with i y .) T and S 1 generate the set of all modular 
transformations. 

The T- and 6> -matrices describes how the set of ground states transform under their respective modular transformations ll48l . 
Generically T, S are elements in a projective representation of SL(2, Z); the double cover of the modular transformations. As 
discussed in the previous section, Tis a diagonal matrix with entries 9 a known as the 'topological spin', the action of rotating 
a quasiparticle type a by 2ir. S ao gives the mutual statistics of braiding a and b around each other. Note that under the modular 
transformations the fluxes transform as T : ($ x , 3L) — > ($> x ,$ y + $> x ), S : ($> x ,$> y ) — > ($> y ,$ x ). As we have discussed, 
$j = is most natural for bosons, but $j = tt is most natural for fermions, so we must return to this subtlety. 



Constraining (or determining) S from J 7 . When t x — 0, we can use S to relate the two flux matrices 



J-y — SJ- X S 1 . 

While this alone cannot be used to solve for S, there are additional constraints, 



(66) 



Sta 



da 

V 



Sab — Sba, 



(67) 



where d a is the quantum dimension for the quasiparticle a, with T> = \/J2 a being the total quantum dimension. -f a is 
the topological entanglement entropy for the MES |3„) defined in the main text. For certain phases such as the Moore- 
Read state and the v = 2/5 Jain state the modular 5-matrix may be determined from these constraints alone. Solving for 
S in the MES basis essentially amounts to diagonalizing T x . 

• Flux sectors and modular transformations. A subtlety in the computing T and S is the interplay of modular transformations 
with boundary conditions (& x , & y ). Since T. S change the fluxes, we need to instead consider a larger Hilbert space for 
which the boundary conditions may take on four possible combinations: ($ x , Q y ) G {(0, 0), (0, tt), (tt, 0), (tt, tt)}, which 
as a shorthand we refer to as PP, PA, AP, AA, respectively. Each of these sectors consist of m linear independent ground 
states, for a total of 4m-dimensional ground state manifold. The PP sector is closed under the action of T and S, but the 
other three sectors mixes under modular transformations ll60l . 

We write T and S as block matrices, where each block is an m x m matrix describing transitions between sectors. The 
order of the four columns/rows are PP, PA, AP, AA. 



T 



' J-pp 










j-PA 












T+ 






T- 





s 



- S PP 












s+ 






s- 












5 AA . 



(68) 



In the minimal entangled basis, each T-submatrix are still diagonal. In the AP and AA sectors ($ K = tt), the formula 
Eq. ( |oT| squared gives the product T + T _ , as two Dehn twists are required to come back to the same wavefunction. In 
other words, Eq. ( |oT| will only give the topological spin h modulo |. 



